Maximal couplings in VT— symmetric 
chain-models with the real spectrum of energies 



Miloslav Znojil 

Nuclear Physics Institute ASCR, 
250 68 Rez, Czech Republic 
e-mail: znojil@ujf.cas.cz 



Abstract 

The domain T> of all the coupling strengths compatible with the reality of the energies 
is studied for a family of non-Hermitian N by N matrix Hamiltonians H^ N ' with 
tridiagonal and VT— symmetric structure. At all dimensions N, the coordinates are 
found of the extremal points at which the boundary hypersurface dT> touches the 
circumscribed sphere (for odd N = 2M + 1) or ellipsoid (for even iV = 2K). 

1 Introduction 

1.1 Non-Hermitian chain models 

In many quantum systems (typically, in nuclear and condensed matter physics), the 
observed spectra can be fitted by the equidistant harmonic-oscillator energies (i.e., 
by E^ H °^ = 2n + 1 in suitable units). An improvement of this fit can be based on a 
perturbatively mediated transition, say, to the popular nearest-neighbor-interaction 
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model with an infinite-dimensional "chain-model" tridiagonal Hamiltonian 

1 a 
b 3 ax 
bi 5 a 2 



(1) 



For the real coupling strengths a n , b m the manifest asymmetry of our Hamiltonian 
i^( Ar ) (with infinite as well as finite matrix dimension N) need not necessarily con- 
tradict the postulates of Quantum Mechanics. This may be illustrated on the so 
called Swanson's model with N = oo [1] or, more easily, on the simplest truncated 
two-dimensional special case of eq. (1), 



H (2) 



Both its eigenvalues E± — 2 ± y/l + ab remain real (i.e., in principle, "observable") 
whenever ab > —1. Inside this "domain of physical acceptability", i.e., for 




(a, 6) G P {2) = {(x,y) | x,y e R, xy > -1 } 



(2) 



these energies are also non-degenerate. This gives the technically welcome guarantee 
that can be diagonalized in a bi-orthogonal basis formed by the two respective 
sets of the right and left eigenvectors | ±) and | ±)) such that 



H^\±) =E±\±), ((±\H^ = ((±| E±. 



(3) 



The diagonalizability is lost on the boundary dV^ [where the basis of eq. (3) becomes 
incomplete] and the reality of the energies is lost everywhere in the open complement 
of©( 2 ). 

In a way discussed thoroughly in our recent letter [2] the two-dimensional model 
(3) proves particularly useful for an elementary explicit illustration of one of the "key 
tricks" which re-assigns the necessary Hermiticity to the similar operators. The goal 
is being achieved by a suitable redefinition of the metric 9 and, hence, of the scalar 
product, 

|V) ©10} = MO|0), = t >O. (4) 

During the last few years, such a recipe has been revealed and/or implemented by 
several independent groups of authors sampled here in refs. [3, 4, 5, 6]. 
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1.2 Construction of the metrics G for a given Hamiltonian 

One should re-emphasize that in general, a redefinition of the metric in Hilbert 
space is fully compatible with the postulates of Quantum Mechanics, provided only 
that the reality of the spectrum is guaranteed. From time to time, the efficiency 
of the trick is being confirmed in various less standard applications of Quantum 
Theory [7, 8]. 

In the notation of eq. (3) the essence of the trick derives from the observation that 
for many manifestly non-Hermitian Hamiltonians H ^ W with real spectra one can 
follow the two-dimensional guidance and construct the two families of the left eigen- 
vectors \n)) and of the right eigenvectors \m) of a given H . They form a biorthogonal 
basis in Hilbert space. In the next step one easily verifies, in all the finite-dimensional 
cases at least, that the operator defined by the spectral-representation-like formal 
expansion 

© = E l n » s n ({n\ (5) 

n 

satisfies the linear operator equation 

& = & H . (6) 

In the final step of the argument one restricts all the parameters s n to the real 
and positive numbers and concludes that the properties of the resulting operator 6 
qualify it for a metric-operator interpretation as discussed in the review paper [3]. 
This means that the "correct" inner product is ambiguous as it may be defined by 
any prescription (5). Its choice in fact fixes our selection of an explicit representation 
of the Hilbert space of states and, hence, "the physics" . 

There exist several remarkable differences between the unique, "standard" choice 
of 6 = I and all the "nonstandard" 6 ^ I in eq. (4). For this reason, usually, 
the Hermiticity condition (6) with 6 ^ I is being re-named to "quasi-Hermiticity" 
[3, 5]. One of the most characteristic consequences of the quasi-Hermiticity of a 
Hamiltonian H lies in the necessity of a specification of the domain T> of parameters 
where the spectrum of energies remains real. 
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2 VT— symmetric models 



2.1 Modified harmonic oscillators 

Under the assumption H ^ W some of the eigenvalues become complex whenever 
we leave the quasi-Hermiticity domain T> of parameters in H. In the context of one- 
dimensional differential Schrodinger operators the problem has been made popular 
by Bender et al [4] who studied the generalized Bessis' oscillators 

H^ GB \v) = -^- 2 +g{x)x\ g(x) = (ix)», v e M (7) 

and conjectured that all the bound-state energies remain real iff v > 0, i.e., inside the 
very large domain V^ GB ^ = (0, oo) of the exponents v. Rigorously, this conjecture has 
only been proved three years later [9]. One should notice that the difficulty of this 
proof is in a sharp contrast with the elementary character of the above-mentioned 
construction of related to the finite-dimensional H^. 

Our present paper is inspired by the question of feasibility of the constructions 
of the quasi-Hermiticity domains T>( N ^ for matrices at the higher dimensions N > 2. 
Predecessors of such a project can be seen not only in the exhaustive analyses of 
virtually all the two-dimensional cases [10] but also in our recent note [11] where 
we reported the feasibility of a complete and non-numerical reconstruction of the 
domain T>^ for certain special VT— symmetric three by three toy Hamiltonians. 

We shall address here the natural question of the specification of D for the matrix 
family of the perturbed harmonic oscillator Hamiltonians (1) restricted by an addi- 
tional requirement of their VT — symmetry. We believe that such band-matrix models 
are really exceptional. One of our reasons originates from the observation that in the 
most elementary differential-operator representation of H^ HO \ all the wave functions 
ip n (x) pertaining to the above-listed energies E^°^ = 2n + 1 are endowed with an 
additional, parity quantum number, V ip n {x) = ip n (—x) = (—l) n i[) n (x) . This is a 
consequence of the commutativity V H^ H °^ = H^ HO "> V which is manifestly broken 
in all the perturbed matrix models H^ N \ 

The VT— symmetry requirement VT = VT is quite natural to impose, 
especially because the operator T can be treated as a mere transposition. In addition, 
it is easy to imagine that in the given basis the operator V is represented by the 
diagonal matrix with elements V nn = (— l) n . As a net consequence, the requirement 
of the VT — symmetry degenerates to the elementary rule a n = —b n at all subscripts 
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n in eq. (1), 



1 











-a 



3 



H (N) = 



5 







(8) 



&N-2 







-ajv-2 2JV-1 



These are the models which we are going to analyze. 

2.2 An additional "up-down" symmetrization 
2.2.1 A generic attraction of the levels. 

After a few numerical experiments with eq. (8) one reveals a comparatively robust 
survival of the reality of the spectrum in perturbative regime. The phenomenon can 
be understood as one of the mathematically most interesting consequences of the 
"sufficient separation" of the matrix elements 1,3,... along the main diagonal [12]. 

In contrast, even a strict observation of the equidistance of the elements on the 
main diagonal need not necessarily be of any help in a deeply non-perturbative 
regime. This danger is well known and the monograph [13] can be consulted for an 
extremely persuasive illustration of the emergence of unexpected difficulties even in 
a symmetric nonperturbative version of our example (8) with an apparently innocent 
choice of the dimension N = 20 and with an apparently "not too nonperturbative" 
constant diagonal where a — a± — . . . — ai 8 = a. 

In order to avoid similar complications in our present VT— symmetric models we 
may try to assume, in the first step, that just a single coupling a = a k becomes large. 
In such a case, solely the two neighboring energies become involved and perceivably 
modified. Generically, in a way controlled by the mere two-dimensional submatrix 



so that the energy values become "attracted" by each other in proportion to \a k \ at 
any k < N — 1, 



mil of H (JV) we have 




E k = 2k + 2 - Vl - a 2 , 



Ek+i 



2k + 2 + Vl-a 2 . 
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The mechanism of this effect is virtually independent of the rest of the spectrum 
(which may be considered pre-diagonalized) so that we may always expect that some 
energies get complexified whenever the couplings become sufficiently strong. 

This means that, intuitively, we may always visualize the coupling-dependence 
of the energies as their mutual attraction. In this sense we are able to guess that 
the levels E no in the middle of the matrix (i.e., such that n « n N/2) will 
be "maximally protected" against the complexification due to their multiple and 
balanced "up" and "down" attraction by all the other levels. 

Such a balance may be quantitatively (though not qualitatively) violated by the 
differences in the absolute values of the pairs of couplings a no+ k and a no _£ at all the 
allowed index-shifts k. For this reason we shall restrict our present attention to the 
special class of the Hamiltonian matrices (8) which are, in this sense, symmetrized 
and have 



dN-2-j j = 0, 1, . . . ,j r , 



entier[N/2] 



(9) 



This means that everywhere in what follows we shall reduce the class of the (N — 
1)— parametric chain models (8) to its "up-down-symmetrized" subset 



H (N) 



-a 






a 
5-3 







5 — 5 ai 















a 

2N + 1 



(10) 



Obviously, the specific choice of the global shift 5 = N of the origin of the energy 
scale makes also the main diagonal of the whole matrix "up-down" symmetric. Still, 
it is slightly unpleasant that at the strictly integer half-dimensions K = N/2, the 
"last" free parameter aj max enters our matrix "anomalously", i.e., just twice. 

This means that the parity of N introduces a fairly nontrivial difference between the 
corresponding up-down-symmetric models (10). 



2.2.2 Even dimensions N = 2K, K = 1,2,.... 

As long as we intend to analyze the secular determinants of our matrices [13], it makes 
sense to simplify our notation and, in particular, to get rid of the subscripts and 
abbreviate a = aj max , b = aj max -i and so on, up to the last element a abbreviated, 
whenever needed, by the last letter z. In this notation, obviously, the symbol z = a 
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coincides with a = ct,j max = a at K = 1 in 



or with b at K = 2 in 





1 


a 








— a 


-1 
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b 










-b 


1 


a 










—a 


-1 


b 










-6 


-3 



etc. We see that the general matrix (10) with iV = 2K may be easily understood as 
partitioned into four K— dimensional submatrices, 



2K - 1 z 








—z ' ■ ■ 








-.3 


b 







: ■■■ -b 


1 


a 


... 


... 


—a 


-1 


b 







-6 


-3 '•• 








z 








-z 1-2K 



The simplest illustrative example has already been shortly discussed above 
(cf. also [2]). In the general case the secular polynomial det^H^ 2 ^ — E) will be a 
polynomial of the K— th degree in s = E 2 and it will only depend on the squares of 
the couplings a 2 max = a 2 = A, a 2 jmax ^i = b 2 = B , . . ., = z 2 = Z . 
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2.2.3 Odd dimensions N — 2M + 1, M — 1,2, 

Whenever the dimension of our band-matrix Hamiltonian with equidistant 

matrix elements on its main diagonal is odd, N = 2M + 1, we have 



H (2M+l) 



2M 


z 

















— z 
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—a 





a 

















—a 


-2 























z 

















— z 


-2M 



Here the central matrix element vanishes and the VT— symmetric coupling is me- 
diated again by the M real matrix elements a,b, . . . ,z. Omitting the overall factor 
E we may reduce the secular polynomial &et(H ( 2M+1 ) — E) to a polynomial of the 
M— th degree in s = E 2 . It will again depend on the squares of the couplings only. 



3 Hamiltonians of the even dimensions TV = 2K 

In the two-dimensional case with K — 1 the whole discussion remains entirely ele- 
mentary (see above) and one can conclude that there exist precisely two points of 
the boundary dV^ (called "exceptional points" in the literature [14, 15]) which are 
defined by the elementary rule a± P ^ = ±1, i.e., by the single root A^ EP ^ = 1 of the 
single energy- degeneracy condition. 



3.1 Four by four model, K = 2 

For the four by four Hamiltonian the standard definition of the spectrum 



det 



3-E 6 

-b 1-E a 

-a -1-E b 

-b -3-E 



i.e., the quadratic secular equation for s = E 2 , 

s 2 + (-10 + 2 b 2 + a 2 ) s + 9 + 6 b 2 - 9 a 2 + b 4 = 



8 



can easily be solved in closed form, 

s = s± = 5 - b 2 - 1/2 a 2 ± 1/2 V64 - 64 b 2 + 16 a 2 + 4 6 2 a 2 + a 4 . (11) 

These formulae may be read as an implicit definition of i.e., of the reality 

domain of the energies or, equivalently [5], of the quasi- Hermiticity domain of the 
Hamiltonian of our K = 2 chain model. 

For a more explicit construction of we can make use of the up-down sym- 
metry (9) and imagine that during the initial perturbative mutual attraction of the 
neighboring levels one can only guarantee the growth of the ground-state minimum 
E = E_ j+ = —^fs+ and the decrease of the top-state maximum E 3 = E +:+ = 

+v /5 +- 

Beyond perturbative domain, at certain "exceptional-point" combinations (a, b) = 
(a,bY EP+ ^ of the sufficiently large strengths a and b, the latter two extreme energy 
levels will ultimately coincide (and, immediately afterwards, complexify) in a way 
discussed in paragraph 2.2.1 above, E^^ = E^P = 0. At another set of the EP 
coupling doublets (a, b) = (a,b)^ EP ~\ both the two "internal" levels may also coin- 
cide as well, E±_ = ^^/sZ = 0. In this way, the complete boundary dT>^ of the 
quasi-Hermiticity domain is a curve in the a — b plane formed by the "weaker" dou- 
blets of the EP-strenghs (a,b)( EP± \ The shape of such a boundary can be deduced 
from eq. (11) (cf. Figure 1). 

In a magnified detail, Figure 2 demonstrates that the graphical representation 
ceases to be reliable in the fairly large vicinity of the common maximum of the 
sizes of the allowed couplings a and b. Fortunately, near any such a "extremely 
exceptional" point (a, b) — { (±V A( EEP \ ±V B ( - EEP ^} of the a - b plane, the details 
of the shape of the boundary dT>^ can be described by the purely analytic means. 
An extension of the latter observation to all the dimensions will become, after all, 
a core of our present message. 

Let us explain the method for = 2K at any K. In the first step one realizes 
that s <yEEF * > = due to the up-down symmetry. As long as this must be the only 
root (i.e., a maximally degenerate root) of the polynomial secular equation 

s K + Pk-i(A, B,...) /- 1 + P K -2(A, B,...) s K - 2 + . . . = (12) 

it can only exist if the K values A <yEEP \ B ( - EEP \ ... of the EEP coupling strengths 
satisfy the nonlinear set of the following K necessary conditions, 

p K ^ (a( eep \b( eep \...)=o, 
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Pk-2 (A& EP \BV EEP \..) = Q, 

(13) 

P (A^ EEP \B^ EEP \...)=0. 
At K = 2 the latter set of polynomial equations reads 

,4 + 25 = 10, (3 + B) 2 = 9A 

and an elimination of A leads to a quadratic equation for B + 3 giving a spurious 
solution A = 64 and I? = —27 (which would imply an imaginary coupling b) and the 
unique correct solution A^ EEF " > = 4 and B <yEEP " > = 3. 

3.2 Six by six model, K = 3. 

In a full parallel with the preceding subsection, secular equation 



det 
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in its polynomial form (12), 

s 3 + (2 b 2 - 35 + 2 c 2 + a 2 ) s 2 + 

+ (6 4 + 2 cV - 44 b 2 + 28 c 2 - 34 a 2 + c 4 + 259 + 2 b 2 c 2 ) s+ 

+a 2 c A - 10 b 2 c 2 + 30 c 2 a 2 + 225 a 2 - 30 c 2 - c 4 - 25 6 4 - 225 - 150 b 2 = 

remains solvable in closed form. As long as our present attention is concentrated 
on the EEP extremes, we shall skip the details of the complete description of the 
hedgehog-shaped surface dT>^ in the full three-parametric space and note only that 
this shape must be all contained within the ellipsoid with the boundary described 
by the first constraint of eq. (13), a 2 + 2 b 2 + 2 c 2 = 35. 

At K = 3 the full solution of the triplet of eqs. (13) ceases to be easy but it still 
remains feasible. Besides the unique and acceptable correct solution 

Al EEP ) = 9, B( EEP ^ = 8, C( EEp 1=5, K = 3, (14) 
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one obtains another set of the alternative solutions generated, after the patient elim- 
ination of A and B, in terms of roots of a final "effective" polynomial in single 
variable C, 

416 C 4 + 20909 C 3 + 22505 C 2 + 28734375 C- 48828125 = 0. K = 3. (15) 

Out of its two real roots, C_ = -65.80360706 and C+ = 1.693394621, the former 
one is manifestly spurious giving the imaginary coupling c. For the latter root we 
have to recall the corresponding condition 

22156250 B + + 2912 C+ 3 + 1446363 C + 2 + 820546875 + 9654410 C + = 

to see that the coupling b = \[B~l is imaginary and should be rejected as spurious 
as well. 

3.3 Eight by eight model, K = 4. 

Out of the four EEP constraints (13) at K = 4 the first equation P^a 2 , b 2 , c 2 , d 2 ) = 
defines the surface of an ellipsoid or, after the change of variables a — > A = a 2 etc, a 
planar side of a simplex, 

A + 2B + 2C + 2D = 84. 

By construction, the domain is circumscribed by this ellipsoid or simplex. 

Unfortunately, one hardly finds any immediate geometric interpretation of the re- 
maining quadratic, cubic and quartic polynomial equations P2(A, B,C, D) = 0, 
P!(A,B,C,D) = and P (A, B, C, D) = containing 13, 19 and 20 individual 
terms, respectively, and admitting just marginal simplifications, e.g., to the 9-term 
equation 

1974 + (B + C + D) 2 + 2 AD + 2 BD + 2 AC = 83 A + 142 B + 70 C - 50 D 
in the P 2 — case, etc. 

In this setting it comes as a real surprise that the above-derived K = 3 rule (14) 
still finds its unique K = 4 counterpart which, in addition, possesses the closed form 
again, 

A ( EEP ) = W, B (EEP ) = 15i C (EEP 1 = 12, D( EEP 1 = 7, K = A. (16) 
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Its derivation necessitated the use of the fully computer-assisted Groebner-basis tech- 
nique. Just for illustration one may mention the K = 4 form of the final "effective" 
polynomial equation, 

314432 D 17 - 5932158016 D 16 + . . . + 153712881941946532798614648361265167 = 0, 

representing the "next-door neighbor" of the still exactly factorizable eq. (15). 

In a test of the uniqueness of solution (16) one finds out that it possesses seven real 
and positive roots D. Out of them, the following three ones are negative and, hence, 
manifestly spurious, -203.9147095,-156.6667001,-55.49992441. We skipped the 
proof of the spuriosity for the remaining four roots, viz., of 0.4192854385, 5.354156128, 1354.675195 
and 18028.16789 since the related calculations, however straightforward, become un- 
pleasant and clumsy. For example, the values of A are given by the rule a x A = (a 
polynomial in D of 16th degree) where the number of digits in the auxiliary integer 
constant a exceeds one hundred. 

3.4 Arbitrary even dimension N = 2K. 

Even though we did not dare to test the applicability of the Grobner-basis technique 
at K = 5, we were lucky in noticing that the previous results already admitted the 
following extrapolation to any K, 

A (EEP) = R 2^ B (EEP) =K 2_ l2> C (EEP) = R 2 _ ^ -q(EEP) = K ^_^^ ^ 

This is our first main result. The validity of this empirically revealed rule has sub- 
sequently been tested and confirmed by the incomparably simpler direct insertions. 

As a byproduct of these verifications, the general elipsoidal surface form of the 
first item in eq. (13) has been predicted from the data available at K < 4 and re- 
confirmed at several higher K > 4 giving, in terms of the original coupling-strength 
variables of eq. (8) with symmetry (9), 

A+2 (B + C+... + Z) = a? +20? + . . .+2 a 2 = V o\ = (18) 

fc=0 

or, in the form of an immersion of T> in an ellipsoid in K dimensions, 

a i + 2b 2 + ... + 2z 2 = J2 4< \ • (19) 
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These observations are in a complete agreement with the individually evaluated for- 
mulae and carry a geometric interpretation showing that every domain X>( 2K ) (where 
all the energies remain real) is circumscribed by a certain elipsoidal hypersurface. Its 
intersections with the boundary dV^ 2K ^ coincide with the 2 K EEP points with the 
coordinates a {EEP ^ = ±K, 6 (£SP) = ±VK 2 - 1 etc. 

4 Hamiltonians of the odd dimensions N = 2M + 1 

In a one-parameteric three-by-three illustration with M — 1, 



2 a 
— a a 
-a -2 



the determination of the interval of quasi- Hermiticity a £ T>^ = (—y/2, \/2) is trivial 
since the secular equation — E 3 + (4 — 2 a 2 ) E = is exactly solvable. In a remark [11] 
we also studied a "generic" three-dimensional (and three-parametric) matrix model 
where we relaxed both the equidistance assumption concerning the main diagonal 
and our present simplifying "up-down" symmetrization assumption a = a 1 . 



4.1 Five by five model, M = 2. 

A comparatively elementary two-parametric example of our present class of models 
of section 2.2.3 is still encountered at M = 2, 



#(6) 



4 

-6 






b 

2 a 

-a 

-a 








a 



-2 
-b 






b 



Its secular equation gives the central constant energy E = 0. The other two pairs 
of the real or complex conjugate levels E n = —E_ n = y/s with n — 1, 2 are obtained 
from the remaining polynomial equation in the new variable s = E 2 , 



-s 2 + (20 - 2 b 2 - 2 a 2 ) s - 64 - 16 6 2 + 32 a 2 - 6 4 - 2 a 2 6 2 



0. 



(20) 



13 



We should determine the domain T>^ in which all the energies remain real. This 
means that inside the closure of the domain of quasi-Hermiticity both the roots 
of eq. (20) must be non-negative. 

Our task is elementary since the M = 2 eigenvalue problem is solvable in closed 
and compact form, 

E ±1 = ± >fw - a 2 - b 2 - V36 + 12 a 2 + a 4 - 3W , 

E ±2 = ±\/ 10 - a 2 - b 2 + V36 + 12 a 2 + a 4 - 3&V 2 . 

Thus, the results of the method of preceding section may be complemented by direct 
calculations. In terms of the two non-negative quantities A = a 2 > and B = b 2 > 
the reality of the energies will be guaranteed by the triplet of inequalitites. The first 
one reads 10 > A + B and restricts the allowed values of A and B to a simplex. 
The second condition 36 + 12 A + A 2 > 36 B requires that the allowed values of 
B must lie below a growing branch of a parabola B max = B max (A). The third 
condition (8 + B) 2 > (32 — IB) A represents an easily visualized upper bound for 
A < A max = A max (B) where the latter hyperbola-shaped function grows with B in 
all the interval of interest. 

Beyond the above direct proof we may also parallel the considerations of the 
preceding section and imagine that the symmetry of eq. (20) implies that its triple 
root must vanish, s = s^ EEP ^ = 0. This means that in the polynomial eq. (20) both 
the coefficients at the subdominant powers of s must vanish. These two coupled 
conditions degenerate to the single quadratic equation with the unique non-spurious 
solution A ( - EEP ^ = 6 and B ( - EEP ^ = 4. Thus, in a way complementing our above 
direct discussion of the reality of the energies we see that at our EEP point all the 
three above-mentioned inequalities become saturated simultaneously. 

4.2 Seven by seven model, M = 3. 

By the same Grobner-basis method as above we derive the result 

A (EEP) = 12 , B (EEP) = 1Qj C (EEP) = Q^ M = 3> (21) 

It is again unique because one of the two roots C± = 27 ± 9 \p2A of the "first 
alternative" Grobnerian "effective" equation C 2 — 54 C = 972 and both the roots 
-354 ± 60 v 7 ^ of the "second alternative" equation C 2 + 708 C + 2916 = arc 
negative while the only remaining positive root C+ = 68.24318125 gives the negative 
B = 28 - 3 C. 
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4.3 All the (2M + 1) -dimensional models with M > 4. 

At M = 4 we still were able to evaluate the explicit form of the secular equation, 

14745600 - 7372800 A + . . . + (-2 C + 220 - 2 B - 2 A - 2 ff - 2 D) s 4 - s 5 = 
and we also still computed the M = 4 EEP solution directly, 

A (eep) = 20 ^ B (eep) = 18 ^ C {EEP ) = IA, £>( BEP ) = 8, M = 4. (22) 

We already gave up the discussion of its uniqueness as overcomplicated. Starting 
from M = 5 this enabled us to change the strategy and to continue the calculations 
by merely confirming the validity of the following general odd-dimensional formula 

a (eep) = M ( M + ^ b (eep) = M ( M _|_ i) _ i . 2 = M(M + 1) - 2, 

c (eep) = M ^ M + !) _ 2 . 3, D( EEP ^ = M(M + 1) - 3 • 4, .... (23) 

This formula is our second main result. 

In order to complete the parallels with the previous section, let us finally recollect 
the universal elipsoidal-surface embedding (18) of the domains £>( 2X ) and emphasize 
that its present odd-dimension analogue is even simpler. Indeed, returning once more 
to all the M < 4 calculations of this section we arrive at the extrapolation formula 

2.U :i • 3.1/ 2 • M 

A + B + C + D + ... + Z = -— — (24) 

o 

the validity of which is very easily confirmed (and was confirmed) at a number of 
higher integers M > 5. Its alternative arrangement reads 

a 2 + b 2 + ... + z 2 < (25) 

showing that every quasi-Hermiticity domain 2}( 2M+1 ) is circumscribed by a certain 
minimal hypersphere, with the mutual intersections lying precisely at the 2 M EEP 
points with the coordinates a^ EEP ^ = ±y / M(M + 1), b^ EEP ^ = ±^M(M + 1) - 2 
etc. 



5 Summary 

We introduced a class of the tridiagonal and up-down symmetrized matrix chain 
models H^ N \ the spectrum of which remains equidistant in the decoupled limit. We 
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believe that beyond their above-mentioned direct connection to physics of harmonic 
oscillators exposed to a small finite-dimensional perturbation, another interesting 
source of their possible future physical applicability could be sought in the equidis- 
tance of spectra of certain manifestly finite- dimensional spin-chain models possessing 
equidistant spectra (cf., for illustration, the Polychronakos' SU(N) model [16] or its 
supersymmetric SU(m|n) generalization [17] etc). 

At any dimension N of our Hamiltonians we determined the coordinates 

of all the EEP (= extreme exceptional point) iV— plets of the matrix elements 
a (EEP) ^(eep) ^ ^ z (eep) ^ choice of which leads to the maximal, iV— fold de- 
generacy of the N— plet of the real energy levels pertaining to the underlying model. 

At N = 2M + 1 the latter EEP values are "maximal" in the sense of the norm 
defined as a square root of the sum of their squares. The same comment applies at 
the even dimensions N = 2K after a slight modification of the norm taking just one 
half of the value of the "central" coupling a 2 in the sum displayed in eq. (18). 

Some of the specific merits of our class of models may be seen 

• in the "user- friendly" tridiagonal structure of its Hamiltonians H^; 

• in the feasibility of an illustrative simulation of all the possible scenarios leading 

to 2k— tuple EP-like degeneracies of the energies (followed by their VT— symmetry- 
related complexifications) at all the eligible multiplicities k < N/2; 

• in the fact that for the latter and similar purposes the models contain precisely 
a necessary and sufficient number of free parameters; 

• last but not least, in an "exact solvability" leading to closed formulae at all the 
dimensions N, for the EEP coordinates at least. 
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Figure captions 

Figure 1. One quarter of the domain T>^ (cf. paragraph 3.1) 
Figure 2. A magnified spike of the domain T>^ 
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